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Abstract 

The generating functional for hard thermal loops in Quantum Chromodynamics is impor- 
tant in setting up a resummed thermal perturbation theory, so that all terms of a given 
order in the coupling constant can be consistently taken into account. It is also the func- 
tional which leads to a gauge invariant description of Debye screening and plasma waves 
in the quark-gluon plasma. We have recently shown that this functional is closely related 
to the eikonal for a Chern-Simons gauge theory. In this paper, this relationship is explored 
and explained in more detail along with some generalizations. 
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1. Introduction 



We have recently shown that the generating functional for hard thermal loops in 
Quantum Chromodynamics (QCD) is very closely related to the Chern-Simons (CS) gauge 
theory ^. In this paper, we shall explore in further detail this remarkable connection 
between QCD at finite temperature and the CS theory. 

Conceptually, the question we are considering is very simple. We consider QCD at 
temperatures well into the deconfinement phase. Thus we have a hot plasma of quarks and 
gluons. As is well known, one expects to have Debye screening in a plasma. If we consider 
an Abelian plasma of positive and negative charges ±e, viz. electrodynamics, screening can 
be understood using the classic argument of Debye. One considers the Poisson equation 
for the electrostatic potential of a test charge, say positive, in the plasma. 



where the right hand side is the charge density in the vicinity of the test charge, n is the 
average number density of particles. The exponentials are the Boltzmann factors giving 
the preferential accumulation of negative and depletion of positive charges in the vicinity 
due to the Coulomb forces, with proper normalization. The approximation (lb), which 
is valid for high temperatures, shows that the solutions have the screened Coulomb form 
iexp(— m£)r) with a Debye screening mass = (ne^/T). For a relativistic plasma, the 
qualitative features of this argument are valid and with n ~ T^, we expect m'j^ ~ e^T^. 
And by calculating the photon propagator in thermal electrodynamics, one can indeed 
obtain this result ^. 

The above argument is presented in terms of potentials and as such, it does not seem 
to be gauge invariant. For the Abelian case, one can easily reformulate the arguments 
using only the gauge invariant electric and magnetic fields. However, in a non-Abelian 
plasma, such as the quark-gluon plasma, it is difficult to avoid the use of gauge potentials 
altogether. Further, since even the notion of the charge of a gluon has to be defined with 
respect to some chosen Abelian direction of the gauge group, it is clear that the simple 




(la) 



(16) 
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argument of equation (1) will have to be modified. The question of interest to us is: how 
do we obtain a gauge invariant description of Debye screening in a non-Abelian plasma (in 
terms of gauge potentials)? More specifically, we need a functional of the gauge potentials, 
T[A], which is generated by the statistical distributions and is effectively a gauge invariant 
mass term for the gauge fields. Screening, of course, is the static part of a more general 
problem, the dynamical part of which is the propagation of plasma waves. Such a r[^] 
will thus be important for the discussion of plasma waves as well. 

Hard thermal loops are closely related to the above. They arise because of the need 
to carry out a partial resummation of perturbation theory in thermal QCD. The need 
for resummation is easily seen using the following argument due to Pisarski ^. Consider 
the elementary polarization diagram of gluons, fig. (la) and its first correction fig.(lb). 
Before the final loop integration over p, the ratio of diagram (lb) to (la) is n(p)/p^. 
Because of Debye screening, the polarization tensor Il{p) has a small p-expansion, Il{p) ~ 
g^T^ + g'^T\p\ +.... ((7 is the quark-gluon or gluon-gluon coupling constant.) We thus see 
that for the small p-regime of integration, i.e. p < g^T^, the naively higher order diagram 
(lb) is comparable to the lower order term (la). Therefore, to be consistent to a given 
order in g, one must sum over (la), (lb), and a series of further insertions of n(p). This 
resummation leads to new effective propagators, and also new effective vertices in general. 
The resummation can be summarized by saying that the new propagators and vertices 
arise from an action 



where we add a term r[A] to the standard Yang-Mills action for gluons. 

Going back to figures (1), we see that the regime of interest is when the momentum 
external to Hip), viz. p is of the order gT or less. This holds in general and the momenta 
carried by the potentials in r[^] can be taken to be ~^T. General power counting argu- 
ments by Pisarski and Braaten show that r[^] is the sum of 'hard thermal loops' ^. These 
are one-loop diagrams where the external momenta are ^gT and the loop momentum is of 
the order of T, i.e. relatively hard. An explicit form of ^[A], the generating functional of 
these hard thermal loops, is necessary for the proper set up of thermal perturbation theory. 
Although of seemingly different origin, r[^] so defined will be the same as the functional 




(2) 
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describing Debye screening and plasma waves. We shall see that it is also closely related 
to the Chern-Simons gauge theory. 

The Chern-Simons (CS) action, it may be recalled, made its appearance in physics 
literature over a decade ago as a mass term for gauge fields in three dimensions ^. Studies 
since then have revealed a number of fascinating features of this action. The Abelian 
version can be used for spin transmutation, converting spin-zero bosons into anyons ^. The 
correlators of Wilson lines in a pure CS theory arc related to the polynomial invariants 
of knot theory ^. Pure CS theory is also related to conformal field theory and Wess- 
Zumino-Witten (WZW) models in two dimensions Chern-Simons-Higgs theories have 
an intriguing set of vortex solutions ^. Finally actions related to the CS action can be 
used for a Lagrangian description of selfdual gauge fields and integrable systems Our 
observation gives another context, viz. the very physical context of the quark-gluon plasma, 
for the understanding of which the CS theory is very useful. 

This paper is organized as follows. We begin with the explicit evaluation of some hard 
thermal loops. The two- and three-point functions are evaluated. These calculations have 
some overlap with the similar calculations of Frenkel and Taylor However we organize 
and orient the calculations towards a more efficient realization of the CS connection. In 
section 3, generalizations to n-point functions using the gauge invariance of r[^] are con- 
sidered. The strategy of using gauge invariance in this way is, to some extent, parallel to 
the work of Taylor and Wong In section 4, we discuss the pure CS theory highlighting 
the eikonal and other features of relevance to the plasma problem. In section 5, we obtain 
r[A] in terms of the eikonal of the CS theory. In fact, equation (86) is the central result 
of our analysis. Plasma waves are briefiy discussed in section 6; non- Abelian plasmons, we 
argue, must be understood as the propagating solutions of the effective action (88). 

2. Evaluation of Hard Thermal Loops 

We begin with the explicit evaluation of some of the hard thermal loops. We shall 
consider one-loop quark graphs with two and three external gluon lines; analysis of these 
diagrams will suffice to abstract many of the features that generalize to the n-point func- 
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tions. As mentioned in the introduction, there are many related but different formahsms 
for deahng with thermal corrections in a field theory We shall use Minkowski space 
propagators with thermal averages for products of creation and annihilation operators. 
The calculation is conceptually the simplest in this formalism. The relevant part of the 
Lagrangian for the quark fields q, q is 

C = qi^-{d + A)q (3) 

where = —it"^A^ is the Lie algebra valued gluon vector potential, are hermitian 
matrices corresponding to the generators of the Lie algebra in the representation to which 
the quarks belong. We shall not explicitly display the quark-gluon coupling constant ^, as 
it is easily recovered at any stage by A — > gA. The one-loop quark graphs are given by the 
effective action 

r = -z Tt log (1 + 5 7 • ^) (4) 

where 

S{x,y) = {T q{x)q{y)) (5) 

is the quark propagator. Tr in (4) includes the functional trace as usual. The two-gluon 
and three-gluon terms in T are given by 

r(2) ='-j dt'xd^y TV [ 7 • A{x)S{x, y)^ ■ A{y)S{y, x)] (6a) 

r('^ = -^ / d^^d^'yd^^T^il ■ A{^)S{x,y)l ■ A{y)S{y,z)^ ■ A{z)S{z,x)] (66) 

where Tr denotes the trace over the color matrices t"' and the Dirac matrices. The quark 
propagator is given by 



where p° = p = p' = {p^,—\p\), and Q{x), of course, is the step function. 

Also 

Op = 1 - Up, (3p = Hp. (8) 
The distribution functions np,np are defined by the thermal averages 

(6t«.^6^'*) = npS'-'S^'l^ (9) 

where (a|"''^, a^'''), (&|"''^, bp'^) are the annihilation and creation operators for quarks and 
antiquarks respectively, a, j3 are spin indices; r, s are color indices. For a plasma of zero 
fermion number, we can take 



1 

rip - Hp - ^po/r _^ 



(10) 



For a plasma with a nonzero value of fermion number, there is a chemical potential 
and correspondingly n^, rip are not equal. In using the expression (7) for the propagator in 
(6), we encounter many terms corresponding to different ordering of the time arguments 
x^,y^,z^. The strategy will be to carry out the time-integrations first, introducing con- 
vergence factors 6=*=^^° , 6=*=*^^° etc.,e small and positive, as required. The integrations give 
energy denominators and bring the result to a form where simplification appropriate to a 
hard thermal loop, viz. the loop momentum being hard (~ T) and the external momenta 
being relatively soft (~ (/T), can be implemented easily. 

Simplification of tlie two-point function 

In using (7) for the quark propagator, we find four terms in F*^^) with > and 
four terms with > x^. Writing 



(27r)^ 

and carrying out the time-integrations we get 



(11) 



OipPq 



q° -k° - ie p^ - q° - k° + ie 



6 



PpPq 





+ -k^ -ie 










Mq 




pO 


+ + fco - ze 


+ q^ + + 




f 









ypu 



q^ + k^ - ie - gO + A;^ + ie 



(12) 



where 



T{p, g) = Tr[ 7 • A{k)'y ■ p7 • A(/c')7 • q 



(13) 



and 



(14) 



In (12), p = q + k. Since p° = |g + A;| ^ q'^ + q - k/q'^ for \k\ small compared to |g|, the 
denominators in (12) involve k ■ Q, k ■ Q' and 2q'^ + k ■ Q, 2q^ + k ■ Q' where 



Q = (i, ^), g' = (i, -^) 



(15) 



The ie's in the denominators in (12) can be let go to zero at this stage. Of course, this 
requires that the kinematics be so chosen that k ■ Q, k ■ Q' and g° are not zero. Making 
this choice and using a, P from (8), we find, for the temperature-dependent part of T^'^\ 



r(2) = -1 J dfi{k) J 



d^q 1 1 
(2^ V V 



/ N T{p,q) _ T{p',q') 



p^-q^ + k^ 



I ^ _ ^ r(p,g') T{p',q) 



+ gO + A;0 



(16) 



We see that the result is linear in the distribution functions, a property that is expected 
to be true in general ^. Equation (16) agrees with Frenkel and Taylor . Some of the 
further simplification of (16) can be made along the lines of their calculation. For \k\ small 
compared to the loop momentum we have 



p^ -q^ -k^ ^-k-Q p^ -q^ + k^ ^k-Q' 
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(17) 



T(p, q) ~ 8?°tr(^i • QA2 ■ Q) 

Tip',q')c,8qhr{Ai-Q'A2-Q') (18) 

T{p', q) ~ Tip, q') ~ Vtr(^i • Q' ■ Q + A^ ■ QA^ ■ Q' - 2A^ ■ A2) 

where Ai = A{k), A2 = A{k') and the remaining trace in the expressions for T's, de- 
noted by tr, is over color indices. The difference of distribution functions can also be 
approximated as — riq ~ 'j^Q ' ^- Using these results (16) simplifies to 



d^q 
(2^ 



tr 



f dnAi-QA2-Q dn A^ ■ Q' A2 ■ Q\ ^^ - 
^dq^ k-Q dqO k-Q' >^ 



n + n 



— {Ai ■ Q'A2 ■Q + Ai-QA2-Q'- 2Ai ■ A2) 



(19) 



We have the result 



(20) 



for any function / of Q, or Q' . We can further use 2Q ■ k = k ■ Q' — k ■ Q. Expression (19) 
then simplifies to 



r(2) 



-i J dfiik) I 



d^q 



tr 



n + n 



q' 



— {2Ai ■ QA2 -Q-Ai- QA2 -Q' -Ai- Q'A2 • Q + • A2) 



(21) 



The angular integration in (21) over the directions of q (or Q) help simplify it further 
by virtue of 



dQ{2Ai-QA2-Q-Ai-QA2-Q'-Ai-Q'A2-Q + 2AiA2) = j dVL {2Ai-QA2-Q') (22) 
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Defining 



^+ = ^- = ^ (23) 



we can write (21) as 



r<^' = -i/*w/|f,^i6tr 



(24) 

This expression becomes more transparent when written in coordinate space and fi- 
nally using a Wick rotation to Euclidean space. We define the Green's functions 



G{xi,X2) 



-ip-{xi—X2) ^4 



P 



P-Q (27r)^ 



/ -ip-(xi— X2) ^4 

In terms of the null vectors Q = (1, q/q^), Q' = (1, —q/q^), we can introduce the lightcone 
coordinates {u,v,x'^) as 

, •u = ^— , Q-x =0 (26) 

where Q = q/q^. We then have Q ■ d = du, Q' • d = dy. We shall also introduce a Wick 
rotation to Euclidean coordinates by the correspondence 

2v 2u z 

d^ = Q-d^2d„ dy = Q'-d^2d, (27) 
The Green's functions in (25) are the continuations of the Euclidean functions 

1 6(^\xJ-x^) 



Ge{xi,X2) 



21X1 {zi — Z2) 



This leads to the correspondence 
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k-Q' 1 1 



i — > 



k-Q TT Z12Z21 

k-Q 1 1 

7 — ~ 

A; • Q TT 2;i22;2i 

1 1 1 

< — > 



k - Q 27rz 2:12 
where Zij — Zi — Zj, etc. Equation (24) for F^^^ can then be written as 



(29) 



d^q 1 



2 ' (277)3 2^0 

2^. ^2 



16 tr 



(n + n) J d'^x A+{x)A_{x) 



2 yrf^^i rf^22 -4+(,ri).4+(.r2) _ /" ,2 -^i -^2 ^-(.'•i)-4-(x2) 
-n TT / d X n IT d X 



TT TT 



^12^21 



TT TT 



^12^21 



(30) 



We define the functional I{A+) by the formula 

(-1) 



2 T c^^-^i d^^n tv{A+{l) - - - A+{n)) 
d X — 



n 

d'x 



TT 



TT 



2 T d^^i d^Z2 tr(.4+(l),4+(2)) 



TT TT 



^12^21 



^12^23 ■ ■ - Znl 
+ --- 



(31) 



We shall show later that I{A^) is related to the eikonal for a Chern-Simons theory. Equa- 
tion (30) for r"^^) can finally be written as 



r(2) = 



d^q 1 



(27r)3 2q' 



K^'M+,A. 



(32a) 



where 



K[A+,A_] = -16 



(n + n) 



d'^x tr [A+{x)A-{x)) + n iTTl{A+) + n z7r/(^_) 



(326) 



A_] in (32a) denotes terms in K which are quadratic in A. I is obtained from / 

by z z. 

Although we have introduced different distribution functions n, n for quarks and an- 
tiquarks, it is only n + n which is relevant at high temperatures. We see that, by virtue of 
J dfl I{A+) = J dQ I{A-), we can write (32b) as 
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K[A+,A_] = -16 



n + n 



J dt^x ir{A+{x)A_{x)) + mI{A+) + mi{A_) 



(33) 



The integral over the magnitude of q in (32) can be easily carried out. With zero 
chemical potential, 

-T2 



r(2) = -1- j dVt j d^x tr(A+ A_) + i7r/(2)(A+) + mi^^\A_) 



(34) 



Simplification of the three-point function 

When expression (7) for the propagator is used in (6b) for the three-point function 
and the time-integrations are carried out we get a number of terms with different types 
of energy-denominators. The dominant terms will have differences of energies p^,g^,r^ 
corresponding to the three propagators so that for a small A;, we get soft denominators. The 
most important terms will have differences of momenta in both the energy denominators. 
There are six orderings of the time labels y^, and for each ordering we get two terms 
with both denominators soft; one of these terms will involve quark distribution factors 
a, j3 and the other involves the antiquark distribution factors a, ^. By looking at the 
exponentials from the propagators with maximal differences of momenta, we can easily 
write down these terms. As in the case of the two-point function we shall neglect the 
ze-factors. The result is then 



r(3) = _ 

3 



diJi{k) r{ki,k2,k3) 



(35) 



T{ki,k2,k3) 



d^q 1 1 1 
(27r)3 V V 2r 

j3pOlqj3j. CXp(3qCXj- 



+ 



(pO-rO + fcO)(gO-rO + fcO) {p^ - q^ + k^){pO - - k^) 

rp( I I /NT fiptSgdlr — OlpOlql3r 



+ 



+ 



{p^ - q^ - kl){q^ - - kl) {p'^ - q^ - k^){p^ - + k^) 
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]}. (36) 



Using the expressions for a, P,a, P in terms of n, n we can further simphfy these terms as 



+ 



- T{p', q', r') [^^o_^o + /,0)[^o_^o_ + 



Uq n 



In these equations, p — q — k2 and r = g + /cs and 

rfM^ = (2n)Hf'Hk. + + *3) g^HI^ (38) 

The numerators in (36,37) are given by 

T{p, q,r) = Tr [ ^ ■ A^^ ■ pj ■ A^^ ■ q^ ■ ^7 • r] ~ 16?°^ tr(^i • QA2 ■ QA^ ■ Q) (39) 

where Ai = A{ki), etc. 

Using the approximation as in (17) for the energy-denominators, we get 

+ 2 .(^. . Q'^. . . + + J^^]}. (40) 

Further, we can approximate the distribution functions np,nr as 

dn , 1 dn 



npc^nq + ^(-^2 ■Q)^nq+ 2^(^2 ■ Q - k2 ■ Q') 

dfi —* —* 1 dn, 

nr^nq + ^{h ■ Q) ~ + 2^(^3 ■ Q' - ks ■ Q) (41) 

When these expressions are used in (40), terms with no ^ factor is zero by momentum 
conservation. With one partial integration over q^ = \q\, the remaining terms simphfy as 
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By using the correspondence rules (29) we can write (42) and F^^) of (35) as 



r(3) 



-167r 



d^q 1 ,2 r d'^zi ^"^^2 d'^zs 
d X 



(27r)3 2(j0 TT TT TT 

tT{A_{xi)A-{x2)A_{x3)) 



n 



tT{A+{xi)A+{x2)A+{x3)) 
^12^23^31 



n- 



^12^23^31 



= / (0^ 2? ^-''^ [ " + ^ z7r/(^)(^-)] , 



(43) 



where I^^^ denotes the terms of cubic order in the A's in the expansion (31) of I{A). Again 
by virtue of J dfl I{A^) = J dfl I (A-), the result depends only on n + n. We can then 
combine equations (32) for the two-point function and (43) for the three-point function as 



r = 



d^q 1 



K[A+,A. 



(44) 



(27r)3 2qO 

where K is given by equation (33) and we retain terms up to the cubic order in A^s. 

We shall now rewrite this result in a slightly different way, with a view to generaliza- 
tions. Using the identity 



J dn ti{A+A_) = -\ 2t^A%AI - J dn {AlAD 



(45) 



and J dO, I{A_) — J dO, I{A^), we can write 



qdq 



{n + n) 



J d^x 27: AlAl - J dn d^x {AlAD - 47ri/(A+) 



(46) 



Notice that F, except for the Aq^q term, depends on A^ only through the combination 
A+ = A-Q/2. (Hereg= |gl.) 

It is possible to interpret the various terms in equation (37) in terms of certain forward 
scattering amplitudes. If we use an interaction Hamiltonian Hint = 7 • ^7 • P and calculate 
the forward scattering amplitude, using standard time-dependent perturbation theory, for 
a quark of momentum p, we get the first term in (37). The sum of this amplitude over 
spin and colors gives the traces. Integration over all momenta, distributed according to np, 
gives the thermal loop contribution. The other terms in (37) can be similarly interpreted in 
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terms of permutations of the external gluons. (This interpretation generahzes to four-point 
and higher functions as well.) 



3. Generalizations 



We shall now consider some of the general features of the calculations presented so far. 
In generalizing to four-point and higher functions, it is easy to see which kind of terms will 
dominate. We need the maximal number of soft denominators. This will come from the 
time-orderings which give maximal number of differences of momenta in the exponential for 
the propagators. In an n-point function, we can get {n — 1) such terms. The denominators 
simplify to products of /c • Q's (and k ■ Q'^s for the antiquark distributions). The general 
result has the structure 



r = — 

n 



+ cyclic 



diJ,{ki, - ■ ■ ,kn) 



(47) 



k2-Q--- k(^n) ■ Q 

and a similar term for the antiquark distributions. (We have renamed A;'s compared to 
(40).) The numerator, in the kinematic regime of interest, viz. A;'s small compared to p, 
becomes Ai ■ QA2 ■ Q ■ ■ ■ ■ Q. We can thus expect the final result to depend on 

A's only through A ■ Q for quarks (and A ■ Q' for the antiquarks). Since J dQ f{A-Q) = 
J f{A- Q'), in fact we can take all higher terms to be a function of A-Q. Further, since 

1 



E 

cyclic 



k2-Qks-Q ■ ■ ■ kn-Q 







(48) 



by conservation of energy and momentum, the nonzero contributions involve ^ and hence 



the quantity is proportional to / p:^^ or T"^ ■ We may then write the result for higher 
order terms in the form 



higher 



r d^i 
~J W 



d^q 1 (n + n) 



r)3 2^0 2 

Combining this with (46), we see that we can write F as 



W{A-Q). 



(49) 



(2 



qdq 



(n + n) 



d'^x 2ttA^A^ + J dn W{A-Q) 



(50) 
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The key observation here is that, including all higher powers of W depends on 
only through A - Q. 

Consider now the contribution of gluon loops. The diagram (4) will generate an Ag- 
term, as in the case of quark loops. For the higher functions, the dominant terms will come 
from the three-gluon vertices. The basic vertex has the form 



where p will be taken as the loop momenta and A denotes the external gluon of momentum 
k. The color structure is given hy A = —iT^^A"' = — /^^A", where /^^ are the structure 
constants of the gauge group. We then see that gluon loops contribute the same way as 
quark loops except that the color trace is over the T"'s which are matrices in the adjoint 
representation. Now, from the structure of the higher terms in A- Q we see that, although 
we have factors like tr(T"iT"2 ■■■T""), all the matrix products can be reduced by the 
symmetry properties of the spacetime part to commutation rules and hence products of 
structure constants except for one last trace over two T's. Thus the ratio of the gluon loops 
to quark loops is given by Ca/Cf where for any representation i?, tr {t°'t'')ji = CrS"'^. For 
gauge group SU (N) with quarks in the fundamental representation, we get Ca/Cf = 2N. 
The fact that there are only two physical polarizations of the gluon is taken care of by 
assigning distribution functions only to the transverse physical polarizations. Equivalently 
we may subtract a ghost-loop contribution. (Naively, (51) leads to 2" tr S^^, = 2" • 4 for 
an n-th order term. Actually, since there are only two physical polarizations we get 2"^+-^ 
as for fermions.) Further there are no separate n and n contributions for gluons. Putting 
all this together we get for the gluon contribution 



Va,i3 = 2iA ■ p6af3 + {^iA ■ k6af3 + 4:ikaA/3) 



2iA ■ p5ai3 



(51) 




(52) 



where Ug is the gluon distribution 



1 



(53) 



n, 



'9 



e9/^ - 1 ■ 
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Strictly speaking we have not proved that W{A • Q) for gluon foops in (52) is the same 
as W{A- Q) of the quark loops in (50). The vertex (51) only shows that the higher order 
contributions, for gluon loops also, depend only on A-Q. The strategy will be to determine 
W{A ■ Q) by gauge invariance of the F. W is thus determined by the gauge transformation 
properties of AqAq and therefore will be the same for quark loops and for gluon loops. 
With zero chemical potential for the quarks we have 

(n + n) TT^T^ 



/ 



qdq 



2 12 
6 



qdq Ug — — - — (54) 



We can combine (52) and (50) for Np flavors of quarks as 



T = {N+INf) 



rp2 



d^x ^ttA^A^q + J dn W{A ■ Q) 



(55) 



127r 

Our arguments so far have been confined to one-loop diagrams. The crucial result 
which can be proved by a power-counting analysis of diagrams is that only one-loop di- 
agrams give hard thermal contributions. Thus in constructing V only one-loop diagrams 
need be considered ^. The other crucial property of F which helps us to proceed further is 
that F, considered as a functional of A^, is invariant under gauge transformations of A^ 
and further that F is independent of the gauge-fixing used for gluon propagators in the 
internal lines. This property can be established by an analysis of the Ward identities. A 
set of identities for the gauge dependence of the generating functional for one-particle irre- 
ducible or proper vertices can be derived using the BRST symmetry. The thermal power 
counting rules, applied to the diagrammatic expansion of various terms in these identities, 
show that the gauge dependence can affect only terms which are of the order of T or 
less, in a high T-expansion. The leading term which is of order T^, viz. the generator of 
hard thermal loops, is gauge invariant and independent of the gauge-fixing chosen for the 
internal lines We shall not repeat this analysis here but notice that since only one- loop 
diagrams are important for F, the thermal contribution is really classical and hence the 
BRST Ward identities can be expected to reduce to a statement of gauge invariance. 

We know that, apart from the AgAo-term, F can be taken to depend on A ■ Q with 
a final integration over the orientations of Q. This result, viz. (55), holds to all orders in 
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A^. Of course, explicitly, we have calculated W{A • Q) to cubic order in A^. We can now 
use gauge invariance along the lines of reference 12, to determine W{A ■ Q). The condition 
for gauge invariance of T is 

J dn5W ^An J d^x A'^uj'' (56) 

where Aq is the time-derivative of Aq and uj = —it°'U!°' is the parameter of the gauge 
transformation, i.e. SA^ = d^u) + [A^,a;]. Equation (56) can be realized by 



6W ^ J d^x A" ■ Qio". (57) 

One can check that (57) is indeed the way gauge invariance is realized by analysis of 
the diagrams. It is clearly so for the two- and three- point functions from our explicit 
calculations. We now rewrite (57), using 

SW = -J d'.(Q.a^^^+lA.Q,j^])V (58) 



as 

where 

We have used the lightcone coordinates from equation (26). We now introduce the fields 

1 5W 

a,^\A-Q^A+ a,^-f = --—-a, (61) 

and also make the Wick rotation, as in (27). The condition of gauge invariance, equation 
(59), then becomes 

dzaz - dzttz + ttz] = 0. (62) 

If 

CL2. •) ^1"^ thought of as the gauge potentials of another gauge theory, we see that 
equation (62) is the vanishing of the field strength or curvature Fzz- The gauge theory 
whose equations of motion say that the field strengths vanish is the Chern-Simons theory. 
We shall therefore turn to a digression on the Chern-Simons theory, returning to (62) 
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and its solution in section 5. Of course, an understanding of Chern-Simons theory is not 
absolutely essential to solving (62). (For a solution in a nonthermal context, see ref.l6.) 
One can simply solve (62) and regard Chern-Simons theory as an interpretation of the 
mathematical steps along the way. However Chern-Simons theory does illuminate many 
of the nice geometrical properties of the final result and is a worthwhile digression. 

4. Chern-Simons Theory 

The Chern-Simons theory is a gauge theory in two space (and one time) dimensions. 
The action is given by 

S=^ [ d^x et''"^ tr(a^a^a« + Wa^Ua). (63) 

471" JMx[ti,tf] 

Here is the Lie algebra valued gauge potential, = —it°'a^. We shall consider SU{N) 
gauge group in what follows. A; is a constant whose precise value we do not need to 
specify at this stage. We shall consider the spatial manifold to be R^, or C since we 
shall be using complex coordinates z = x + iy, z = x — iy. (Actually, we have sufficient 
regularity conditions at spatial infinity that we may take M to be the Riemann sphere.) 
The equations of motion for the theory are 

i^M- = 0. (64) 

The theory is best analyzed, for our purposes, in the gauge where ao is set to zero. In 
this gauge, the equations of motion (64) tell us that a^, are independent of time, but 
must satisfy the constraint 

Fzz = d^az - dzttz + [a^, a^] = 0. (65) 

This constraint is just the Gauss law of the CS gauge theory. It can be solved for Qj;^ ctS Ct 
function of a^, at least as a power series in a^. The result is 

__v^/'_i\n-i f d'^Zn azjzi, Zi)az{z2, Z2) ■ ■ ■ az{Zn, Zn) , . 

- 2.(-i) J - — (^-_^-,)(^-,_^-,)...(^-^_^) • ^^^^ 

This can be easily checked using dz{j^) = 7r5^'^\z — z'). 
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In the ao = gauge, the action becomes 

5 = ^ y dtd^x tr(a^aoa^). (67) 

This shows that is essentiaUy canonically conjugate to a^- In fact in carrying out a 
variation of S, we find the surface term 0{tf) — 9{ti), where 



ik f 

— / (fx tT{a^Sa^). (68) 



'M 

(We assume a^^a^ to vanish at spatial infinity.) 9 is the canonical one-form of the CS 
theory. (This is so by definition; the canonical one-form in any theory can be defined 
by 9, where the surface term in the variation of the action is ^/ — 9i, the subscripts 
referring to the final and initial data surfaces. This is an old result going back to the 
theory of canonical transformations and Hamilton-Jacobi theory. In more modern times, 
for example, Schwinger's action principle is essentially based on this statement. For some 
recent references, see ref.l7. ) ^ is the analogue of Pidx^ of point-particle mechanics; for 

an action S = J dtdx ^(2^)]) we would find 9 = mXiSx^ = PiSx'^. We can make 

another variation of 9, antisymmetrized with respect to the variation Sa^, denoted by the 
wedge product sign, and write 

ik f 

uj = 59 = — / d'^x iT{5ag A Sa^) 

= I /.W.' A (69) 



where 



, ik f 5(x-x')5"^'i rrrc.^ 

Ua,{x,x) = --[^^^_^,^^^, ) 



and 5^" = (5a|,5a"). uj defined by (69,70) is called the symplectic structure and is of 
course the analogue of dpidx^ of particle mechanics. The inverse of ujab gives the Poisson 
brackets, the commutators being i times the Poisson brackets. For our case we get 



or 



[a^(x),a^(x)] = -^rt(2)(x-xO. (71) 
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One does not have to go through the symplectic structure to arrive at (71). One could 
simply use the fact that, from (67) the canonical momenta are tt = a^, 7f = 0. This is thus 
a constrained system in the Dirac sense and using the theory of constraints one can derive 
(71). The procedure of using ujab{x,x') is quicker. 

In the expression (68) for 9, is independent of a^. We can however express as a 
function of via the constraint (65) or equivalently (66) and functionally integrate $. In 
other words, we define /(a^) by 

(5/ = ^ y (fx ti[a^{a^)6a^]. (72) 

The solution for / is given by 

^ n J TT TT ^12^23 •• -^n-ln^nl 

The quantity / has a rather simple interpretation. For one-dimensional point-particle 
mechanics, 6*, as we mentioned earlier, is given by pdx. p is independent of x to begin 
with, but we can express it as a function of x via a constraint such as of fixed energy, 

2 

e.g. ^ + V{x) = E. Integral of 6* = pdx then gives Hamilton's principal function or the 
eikonal, familiar as the exponent for the WKB wave functions of one-dimensional quantum 
mechanics. We have an analogous situation with (65) expressing as a function of a^. / 
is thus an eikonal of the Chern-Simons theory. 

/ is in fact the Wess-Zumino-Witten action We can write the gauge potential 
as — —dzUU~^ where U is in general not unitary; it is an SL{N, C) matrix for gauge 
group SU{N). Notice that since dg has an inverse by virtue of dzjj^^ = 7r5^'^\z — z'), 
such a U can be constructed for any a^, at least as a power series in a^. / can then be 
written as / = —ikSwzwiU) where 

Swzw{U) = ^ I d^x tT{dzUd,U-^) [ d^x e^'^^^iriU-^d^JJU-^duUU-^dJJ). 

27r Jm 127r Jms 

(74) 

As usual the second term involves an extension of U into a three-dimensional space. We 
take = M X [0,1] with U{z,z,Qi) = 1, U{z,z,l) = U{z,z). The relationship to / is 
easily seen by considering variations of Swzw- Under the variation U e'^U ~ {l + (p)U, 
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we find Sa^ = —Dg(fi = —{dg,^ + [a^, (p\) and 

1 



5Swzw = — d X tT{dz(fiaz). 

TT 



(75) 



Partially integrating and using Fg^ = and = —Dz(f, we find that Swzw obeys (72) 
except for a factor {—ik), thus identifying / = —ikSwzw- 

Another quantity of geometrical significance for the CS theory is the Kahler potential. 
The phase space of the CS theory, viz. the function-space of the potentials a^, is actually 
a Kahler manifold, i.e. it is a complex manifold with a metric or distance function defined 
in terms of a potential. Specifically, 



k f 
271" Jm 



d xd X 



2 „/ 



6a1{x) 5al{x') 



K 



5a%{x)5al{x'). (76a) 



The potential functional K for the metric is the Kahler form. Alternatively we may define 
Khj 



d^xd^x' 



Sa^{x) Sal{x') 



K 



5a%{x) ^5a\{x'). 



We find that 



K 



k 

2^ 



d X a%al + h{az) + /i(a^) 



(766) 



(77) 



where h is an arbitrary function of the argument indicated. In general, there is nothing to 
dictate any preferred choice for /?.; but for the CS theory, because there is the additional 
structure of gauge transformations, we can choose h so as to make K gauge invariant, h 
is then proportional to the eikonal / and the gauge invariant Kahler potential is 



K = -- 



1 



TT 



■I 



k I d^x tj:{azaz) + ml{az) + iT:i{az) 



(78) 



Finally notice that the eikonal / of (72) may be regarded as the expansion in pow- 
ers of Qz of the logarithm of the functional determinant of Dz = dz + az', i.e. / = 
{—ik)logdet Dz = —ik TrlogD;^. The Kahler potential K of (78) is then given by 
— /cTr log(-D2-D2). The expansion of this expression in powers of the potential obviously 
gives / and /. The extra term j ^iv{azaz) is precisely the local counterterm needed to 
give a gauge invariantly regulated meaning to Trlog(Dz£)z) . 



21 



5. The Action for the Hot Quark-Gluon Plasma 



We now return to equations (61,62) for the quark-gluon plasma. We can rewrite (61) 



as 



(79) 



Comparing with (72), we see that, since o^, az obey the constraint (62), the solution 
is related to the eikonal /. The difference here is that az = and depend on all 
four coordinates Xf^, not just z, z. However, there are no derivatives with respect to the 
transverse coordinates x'^ in (62) and hence the solution for in terms of is the same 
as in (66), with depending on x'^ in addition to z, z. The argument of all factors is 
the same for x'^, i.e. 

^_ ^ y^._^xn-i f d'^^n az{zi,zi,x'^)az{z2,Z2,x'^) . . .az{Zn,Zn,X^) 

^ ^ J 7T 7T {z-Zi)zi2...Zn-ln{Zn-z) 

For the eikonal we get the same expression as (73) but with integration over the transverse 
coordinates, i.e. 



n 



TT 



TT 



^12^23 • • • Znl 



(81) 



Since it is not relevant to the present discussion, we have, for the moment, set = 1. We 
then find the solution to (79) as 



W = -4:-niI{az) - J alaldt^x. 



(82) 



From (55), F is given by 



T={N+\Nf) 



127r 



dtx 27rylgylg 



(83) 



We can now use identity (45) in reverse to write 



T = -{N+\Nf)— j dVL 
= {N+\Nf)^ j dVt K{A+,A_) 
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d'^x tr(A+A_) + mI{A+) + mi{A_) 



(84a) 
(846) 



where K is given by (78) with the additional integration over coordinates x'^ transverse to 
Q. If we write = —dgUU~^ and A_ = U"^ ^0^17, we can write (84) in terms of Swzw 

as 

r = -(iV + INf)^Swzw{U^U). (85) 

(Again, a suitable additional integration over the transverse coordinates is understood.) 

From equation (55) or (59) it may seem that our solution is ambiguous up to the 
addition of a purely gauge invariant term. But for hard thermal loops, the additional 
structure that W depends only on = \A ■ Q tells us that the gauge invariant piece 
must obey D^^ = {i. Since is invertible, at least perturbatively, there is no nontrivial 
solution to this equation. Thus (83) or (84) is the unique solution and F so defined must 
indeed be the generator of hard thermal loops. 

In the last section, we noted that K{Aj^, A^) can be considered as Tr\og{DzDz)- 
Since and are the chiral Dirac operators in two dimensions, TAog^D^Dz) is clearly 
the photon mass term of the Schwinger model, for Abelian gauge fields. More generally, 
for non-Abelian fields as well, we can consider TAog{DzDz) as a gauge invariant mass 
term. It is perhaps fitting that the gauge invariant Debye screening mass term in four- 
dimensional QCD is given by suitable integrations of such a two-dimensional mass term 
(with, of course, the additional x^-dependence). 

We close this section with two remarks on F. The CS theory (63) violates parity. 
We see that this parity violation disappears, as indeed it should, by integration over the 
orientations of Q, for the QCD case. Alternatively, expression (84a) is manifestly parity- 
symmetric with Aj^ ^ A_, Q ^ Q' under parity. Secondly, instead of setting k = 1 
and then having a prefactor (A^ -|- ^A^f)^ in equation (84), we could simply choose 
k= {N + ^Nf)^. For WZW actions Swzw{U) for which U has a non-Abelian unitary 
part, any action we use must be an integer times Swzw- This has to do with the fact 
that there are homotopically nontrivial maps from the three-sphere into the space of 
matrices U, characterized by the fact that the third homotopy group lis of the space of 
matrices U is the set of integers In our case, F involves SwzwiU'^U) or Swzw{H) 
where H is hermitian. The space of hermitian matrices has trivial lis and so, as expected, 
there is no argument for quantization of the coefficient, which is (N + |A^f)^ for us. 
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We have given the coeflBcient of K in (84) for a plasma with zero chemical potential. 
For the more general case, we can write F as 



Hi + Hi 



1 d^q 



d'^x tr{A+A-) + iTvI{A+) + iTvI{A. 



(86) 



6. Plasma Waves and Debye Screening 



We have obtained the effective action F which is the generator of hard thermal loops. 
As is standard in other contexts of resummations of perturbation theory, we can develop 
the resummed perturbative expansion of thermal QCD as follows. We introduce a splitting 
of the Yang- Mills action as 

S = So-cr 

So = J +F (87) 

We define propagators and vertices and start off the perturbative expansion using Sq. cF 
will be treated as a 'counterterm', nominally one order higher in the thermal loops than 
5*0. Eventually of course c is taken to be 1, so that we are only achieving a rearrangement 
of terms in the perturbative expansion. (As usual we must have gauge fixing and ghost 
terms. We have also not displayed the quark terms; for the quark terms, see ref.20.) Using 
this procedure one can calculate quantities which require resummations such as the gluon 
decay rate in the plasma ^. 

F also gives us an effective action for low energy gluon fields 

Seff = J -If' +r. (88) 

(Although Seff looks like in (87), the interpretation is different. = 5'o — cF in (87) 
generates the thermal perturbation theory. Momenta for fields in Sq, for example, can be 
very high. Calculations starting with (87) lead to a low energy effective action (88); the 
latter is, of course, useful only for small momenta.) 

Long wavelength and low frequency plasma waves are the classical solutions of the 
effective theory (88). It also includes effects such as the screening of Coulomb fields. 
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These features can be seen by examining the Abehan case or electrodynamics. In this 
case, the terms in T which are cubic or higher order in are zero and in terms of the 
Fourier components of A^, we can write 



Seff = i J A^{-k)M^'^{k)A,{k) 



(89a) 



where 



2rp2 



127r 



J k-Q 



(8%) 



k^M^jj = in accordance with the requirement of gauge invariance. We have restored the 
coupling constant e at this stage. We can split A^ into a gauge dependent part and gauge 
invariant components as 

A^ = k^k{k) + + (90) 



where A shifts under gauge transformations and cu^, /^^ are gauge invariant. We take 

( \A ^0 (3) ( \A 



Po = Pi = ei'-^ax, A =1,2 



(91) 



(Here we are considering fields off-shell and so ^ \k\.) The e^'s form a triad of spatial 
unit vectors which may be taken as 



ki 



i = 1,2,3, 



,(1) - 



= (e 



k. ^^^^ ^(,)^^ k,k 



Ixirp 



U2 U2 

tXirptXi 



ih.2 

f), ^=l,2, 



(92) 



where fc^ = fcf + fc^. Notice that kie\^^ = 0, A = 1, 2. (f) and a\ are the gauge invariant 
degrees of freedom in (90). When the mode decomposition (90) is used in (89) we get 



k^ 



)M' {k)Pj{k) + (j){-k)M'^{k)(j){k) 



(93) 



25 



where 



P k^ 2\k\ 



M\k)^k^ + ^^^{1-\L) (94) 
3 ' 2\k\ 

. = log(^^). (95) 

We have also included an interaction term with a source in (93); i.e. we include / A^^J'^ 
and simplify it using (90). From (93) we see that the interaction between charges in 
the plasma is governed by (M^)~^ which shows the Debye screening with a Debye mass 
tud — \f^^^^- The action (93) can also give free wavelike solutions. The dispersion 
rules for these plasma waves would be = for the transverse waves and — for 
the longitudinal waves 

Non- Abelian plasma waves can be defined in a similar way as propagating solutions 
to the equations of motion given by (88). If one truncates T to the term quadratic in A^, 
these are essentially the same as the Abelian plasma waves with mj^ = {N + \Nf)^-^ 
and {N + \Nf) ^ q in M'^{k) rather than ^^^^ . However such a truncation of F loses 
the full gauge invariance. In the approximation of small gA^, i.e. oscillations of very small 
amplitudes, this lack of gauge invariance may not be very serious. However, in general, 
one must really seek propagating solutions to the equations of motion keeping all of F. 
These will be the genuine plasma oscillations in the non- Abelian case. We do not have any 
such solutions as yet, but some of them might coincide with those found by Kajantie and 
Montonen 
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